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1 .  INTRODUCTION 


The  aim  of  this  paper  i-s— to  emphasize  that  there  exists  a  unified  ap¬ 
proach  for  solving  initial  value  problems  for  equations  in  1,  1+1  (i.e., 
one  spatial  and  one  temporal),  and  2+1  (i.e.,  two  spatial  and  one  tem¬ 
poral)  dimensions.  Furthermore  it  remarks  on  inverse  problems  in  higher 
than  two  spatial  dimensions.  Although  these  inverse  problems  are  not 
related  to  physically  significant  nonlinear  evolution  equations,  they 
are  useful  in  the  context  of  inverse  scattering.  In  this  presentation 
we  emphasize  the  main  ideas.  The  detail  formalisms  can  be  found  in 
the  cited  papers. 

It  turns  out  that  solving  the  initial  value  problem  for  some  equa¬ 
tions  for  q(t),  or  q(x,t),  or  q(x,y,t)  is  equivalent  to  solving  an 
inverse  problem  for  some  related  eigenfunction  T(z;t),  or  T(z;x,t), 
or  f(z;x,y,t).  The  inverse  problem  takes  the  form  of  a  Riemann-Hi lbert 
(RH)  problem  for  equations  in  1  and  1+1,  and  the  form  of  a  nonlocal  RH 
problem  or  of  a  5(DBAR)  problem  for  equations  in  2+1  (a  DBAR  problem 
is  generalization  of  a  RH  problem).  To  define  the  relevant  RH  or  DBAR 

•This  article  consists  of  expanded  material  of  six  lectures  presented 
by  one  of  us  (A.  S.  Fokas)  at  the  Workshop  on  “Solitons, "  Ti ruchi rapa 1 1 i  , 
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problems  one  needs  to  study  the  analyticity  properties  of  ?  with 
respect  to  z.  Furthermore  those  problems  are  uniquely  defined  in  terms 
of  certain  asymptotic  data  of  the  underlying  linear  system  satisfied 
by  f:  Monodroay  data  in  1,  scattering  data  in  1+1  and  some  cases  of 
2+1,  and  inverse  data  in  some  cases  of  2+1.  We  use  the  Pamleve  IV(PIV), 
modified  KdV  (mKdV)  and  the  Davey-Stewartson  (DS)  as  illustrative  exam¬ 
ples  of  equations  in  1,  1+1,  and  2+1  respectively. 


The  above  inverse  problems  can  be  naturally  generalized  to  higher 
than  two  spatial  dimensions.  For  example,  the  generalization  of  the 
inverse  problem  associated  with  the  DS  equation  leads  to  an  inverse 
problem  for  a  matrix  valued  function  f ( z ; x^ , x ) ,  z  c  Cn,  XqC**  * 
x  c  Kn,  n  >  1.  However,  while  the  associated  potential  q(Xg,x)  depends 
on  n+ 1  variables ,  the  inverse  data  T(zR,i1,«ij,...,i«n),  ZR  c  '  *  ,  c  . 

e®,  depends  on  3e-1  variables .  This  has  important  implications:  (a) 

The  inverse  data  must  be  appropriately  constrained.  This  “characteri¬ 
zation"  of  the  inverse  data  is  conceptually  analogous  to  the  charac¬ 
terization  of  the  inverse  scattering  data  in  the  multidimensional 
Schrodinger  equation  (1).  (b)  The  existence  of  "redundant"  scattering 

parameters  can  be  used  to  reduce  the  above  problem  to  one  in  two  spatial 
dimensions.  This  is  in  contrast  to  the  case  of  the  multidimensional 


Schrodinger  equation  where  the  inverse  problem  can  be  solved  in  closed 
form.  (c)  Since  the  inverse  problem  for  T  is  reduced  to  one  in  two 
spatial  dimensions,  it  follows  that,  if  one  allows  T  ,  q  to  depend 
parametrically  on  t,  q(Xg,x,t)  satisfies  an  evolution  equation  reducible 
to  two  spatial  dimensions.  In  particular,  the  N-wave  interaction  equa¬ 
tion  in  n+1  spatial  dimensions  can  always  be  reduced  to  two  spatial 
dimensions.  Thus  a  genuine  three-spatial-dimensional  nonlinear  evolution 
equation,  related  to  an  inverse  problem,  remains  to  be  found.  (It  should 
be  noted  that  several  other  "multidimensional"  problems  can  be  reduced 
to  one  or  two  spatial  dimensions,  see  M.  J.  Ablowitz's  contribution 


in  these  proceedings.) 


We  first  define  the  standard  RH  and  DBAR  problems. 


2.  RH  AND  DBAR  PROBLEMS 


Let  C  be  a  simple,  smooth  closed  contour  dividing  the  complex  z-plane 

into  two  regions  D*  and  0  (the  positive  direction  of  C  will  be  taken  0fl 

as  tnat  for  which  D*  is  on  the  left).  ;Av,'“ 
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A  function  $(z)  defined  in  the  entire  plane,  except  for  points 
on  C,  will  be  called  sectionally  holomorphic  if:  (i)  the  function  $(z) 
is  holomorphic  in  each  of  regions  D+  and  D  except,  perhaps,  at  z  =  ■>; 

(ii)  the  function  ♦(z)  is  sectionally  continuous  with  respect  to  C, 

approaching  the  definite  limiting  values  ♦  *(<;)»  ♦  as  *  approaches 

a  point  ^  on  C  from  D* ,  or  D  ,  respectively.  The  classical  homogeneous 
RH  problem  is  defined  as  follows  (2).  Given  a  contour  C,  and  a  function 
G<c )  which  is  Holder  on  C  and  det  G(C !  ♦  0  on  C,  find  a  sectionally 
holomorphic  function  *(z),  with  finite  degree  at  “  ,  such  that 

♦  ^ ( C  )  =  G(;  )♦'<<; ) ,  on  C,  (2.1) 

where  ♦  *(£)  are  the  boundary  values  of  ♦ ( z )  on  C.  If  G(C )  is  scalar, 
(2.1)  is  solvable  in  closed  form.  If  G(C )  is  a  matrix  valued  function, 
then  (2.1)  is  in  general  solvable  in  terms  of  a  system  of  Fredholm 
integral  equations.  Various  generalizations  of  the  above  RH  problem 
are  possible.  For  example:  (i)  The  contour  C  may  be  replaced  by  a 
union  of  intersecting  contours,  (ii)  G(C)  may  have  simple  discontinui¬ 
ties  at  a  finite  number  of  feints;  in  this  case  one  allows  ♦  f  ■» )  to 
have  integrable  singularities  in  the  neighbourhood  of  these  points. 

(iii)  RH  problems  may  be  considered  in  other  than  Holder  spaces  (e.g. 
[3]):  (iv)  One  may  consider  inhomogeneous  RH  problems  ♦*(£)  ■  G(C)$  (C) 
*  F(c )  on  C. 

The  DBAR  problem  can  be  defined  as  follows-  Given  3$/3z,  find  ♦. 

If  3V2z  *  0  everywhere  except  on  a  curve,  then  the  DBAR  problem  reduces 
to  a  RH  problem  (since  3*/3z  ■  -  ♦  ,  in  a  distribution  sense).  The 

DBAR  problem  can  be  solved  via  the  following  extension  of  Cauchy's 
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It  is  interesting  that  the  first  RH  problem  was  formulated  in  connec 
tion  with  an  inverse  problem  (see  [12]  for  references).  Actually  RH 
problems  are  intimately  related  to  solutions  of  inverse  problems  in 
1*1,  2+1,  and  1  dimensions: 


3.  INVERSE  PROBLEMS  IN  1+1 

We  recall  that  a  necessary  condition  for  a  given  nonlinear  equation 
for  q(x,t)  to  be  solvable  via  1ST  is  that  this  equation  is  the  compati¬ 
bility  condition  of  a  Lax  pair  of  linear  equations.  Let  us  consider 
the  modified  KdV  equation 

2 

q+q-6qq*0,  (3.1) 

^t  ^xxx  ^  ^x 

as  an  illustrative  example.  Equation  (3.1)  is  the  compatibility  condi¬ 
tion  of 

Y  (  2  ;  x  ,  t )  =  l  z  (  J  ,  Y  (  z  ;  x ,  t )  ]  +  QY  (  z  ;  x , t )  ; 


1 
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(3.2a) 


Y  t ( 2 ; x , t )  *  [Uq , Y ( z; x,t ) i  +  QY(z;x,t) 


(3.2b) 
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l>4qz2-2iqxz  +  2q3-qxx  2iq2z 


We  first  note  that  the  above  Lax  pair  is  isospectral,  i.e.,  *  0 . 

Aljo  it  turns  out  that  equation  (3.2a)  is  of  primary  importance;  equa¬ 
tion  (3.2b)  plays  only  an  auxiliary  role.  To  solve  the  initial  value 
problem  for  initial  data  decaying  as  |  x  |  ,  one  first  formulates  an 

inverse  problem  for  Y(z;x,t):  Given  appropriate  scattering  data 
reconstruct  Y. 


■  -  -  -  .  ^ 


By  studying  the  analytic  properties  of  f  with  respect  to  z,  where 
t  satisfies  (3.2a)  one  establishes  the  existence  of  a  j  which  is 
a  sectionally  meromorphic  function  of  z,  with  a  jump  along  the  Re  z 
axis.  This  jump  as  well  as  the  residues  of  the  poles,  are  given  in 
terms  of  appropriate  scattering  data.  Thus  the  inverse  problea  is 
equivalent  to  a  matrix,  regular,  continuous  RH  problem  defined  along 
the  Re  z  axis  and  uniquely  specified  in  terms  of  scattering  data. 


Since  in  the  above  discussion  we  have  only  used  (3.2a),  it  is  evident 
that  one  may  pose  an  inverse  problem  for  any  function  q(x).  However, 
this  result  is  useful  for  solving  the  initial  value  problem  for  q(x,t) 
only  if  q  evolves  in  such  a  way  in  t,  that  the  scattering  data  is  known 
for  al-1  t.  If  f  evolves  in  t  according  to  (3.2b)  (i.e.,  if  q  solves 
(3.1).)  then  it  turns  out  that  the  evolution  of  the  scattering  data  with 
respect  to  t  is  simple.  Hence,  the  above  RH  problem  is  specified  in 
terms  of  initial  scattering  data;  its  solution  yields  f(z;x,t)  and 
then  (3.2a)  gives  q(x,t). 
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We  summarize  the  results  of  [  5  jVconcerning  mKdV  in  the  case  of 
solitonless  potentials. 


Proposition  3.1  (Bounded  eigenfunctions).  A  solution  of  (3.2a)  bounded 
for  all  complex  values  of  z  =  z^  +  iz^  and  tending  to  I  as  z  -»<»  is 
given  by 


T  (  z  ;  x ) 


{  . 


t  ( z ; x ) ,  zr  >0 


f  (  z  ;  x )  ,  z  l  <  0 


where  f(z;x)  satisfy  the  following  integral  equations; 


(3.3) 


f  *  (  z  ;  x )  =  I 


iz(  x-f; )  J 


i»tQ(^)f(z;C) 


r  x 

dc< 

-  [  dC  elz(x^)J(w;  *  1’0)Qft 

J  V 
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where  if  F  and  Y  are  2*2  matrices  then 
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*0F  -  Diag(F1  rF22). 


Proposition  3.2  (Departure  from  Holomorphicity-Scattering  Equation). 

f*,  f  are  holomorphic  functions  of  z  for  >  0,  zI  <  0  respectively. 
The  departure  from  holomorphicity  for  z  ■  zR  is  given  by 

A 

f+(z;xJ  -  f” ( z ; x )  -  Y^e1 2xJ ( I-B- 1 < z ) b( z > )  (3.6) 

where 

B ( z )  +  |  dSe"l2^Jfl+(Qf't'),  b(z)  *  I  ♦  f  dCe'lz<-Jir_  (QT_ ) 

J  -flO  ' 

SO, 

»>(z;xlenW(a',|2)b(zl)  »  Y-(z;x).  (3.7) 


Proposition  3.3  (Inverse  Problem-Reconstruction  of  Q) 
Q(x)  is  obtained  from 


Q(  x )  *  (J,  yj-  J  dz’T(z'  ;x)e12  'xJ(  I-B-1  (  z  '  )b(z’  )  )  1  (3.8) 

where  T(z;x)  solves  the  following  Riemann-Hi lbert  boundary  value  problem 


Y  (  z  ;  x ) 


I  + 


1  f  dz  '  T  (  z  '  ;  x  )  e 

1,1  J— 


eiz ' xJ ( I-B~ 1 ( z ‘ ) b( z ' ) ) 


(3.9) 


(z  -  iO) 


Using  equation  (3.2b)  we  obtain: 

Proposition  3.4  (Evolution  of  Scattering  Data).  The  evolution  of  the 
scattering  data  from  B(z;0),  b(z;0)  is  g.ven  by 

Unt  unt 

B(z;t)  -  e  u  B ( z ;  0 )  ,  b(z;t)  -  e  U  b(z;0). 


Since  B  (resp.  b)  is  a  strictly  upper  (resp.  lower)  triangular  matrix 
the  evolution  of  the  scattering  data  is  given  by 


B  |  2  (  z ;  t )  -  e 


b21 ( z ; 0  )  . 


(3.10) 
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4.  INVERSE  PROBLEMS  IN  2+1 

Let  us  consider  the  Davey-Stewartson  equation  (a  two  dimensional  ana¬ 
logue  of  the  nonlinear  Schrodinger  equation) 


iQt  ♦  T<°2Qxx+ayy>  -  -a2MQ|2Q  *  *0.  tx>:  -  02*yy  *  2Xo 2 £  |  Q | 

X  *  t  1 

as  an  illustrative  example.  A  Lax  pair  for  (4.1)  is  given  by 


'1  O' 

'0 

Q  ' 

■  iz(JT-fJ)  +  qf  +  oJT..,  J  * 

.  q  * 

y 

Tt  =  A3fyy  *  A2fy  *  V  '  ^ ( *3*  A30  >  *  2izA3ty  +  izA2T' 

where  AJ(  A2 ,  A^,  A^0  are  appropriate  matrix  functions  of  Q,  Q 
denotes  complex  conjugate). 
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(4.1) 


(4.2a) 


(4.2b) 

(The  bar 


The  situation  is  conceptually  similar  to  the  case  of  1+1:  To  solve 
the  initial  value  problem  for  q(x,y,t)  one  first  formulates  an  inverse 
problem  for  T(z;x,y,t).  Depending  on  the  value  of  o  there  exist  two 
different  cases  (for  brevity  of  presentation  we  assume  non-existence 
of  poles,  i.e.,  non-existence  of  lumps):  (i)  o  *  1.  There  exists  a 
T  which  is  a  sectionally  holomorphic  function  of  z  and  which  has  a 
jump  along  the  Re  z  axis.  This  jump  is  also  given  in  terms  of  scattering 
data  but  it  depends  on  them  in  a  non-local  way.  Thus  the  inverse  pro¬ 
blem  is  equivalent  to  a  non-local,  matrix  continuous  RH  problem  defined 
along  the  Re  z  axis  and  uniquely  specified  in  terms  of  scattering  data., 
(ii)  o  »  -  i.  There  exists  a  T  which  is  bounded  for  all  complex  z, 
but  which  is  analytic  nowhere  in  the  complex  z  plane.  However,  its 
departure  from  holomorphicity  3T/3z  can  be  expressed  in  terms  of  appro¬ 
priate  inverse  data.  Thus,  now  the  inverse  problem  is  equivalent  to  a 
3  (DBAR)  problem:  Given  3T/3z  reconstruct  f  . 


Using  (4.2b),  again  one  shows  that  the  inverse  scattering  and  the 
inverse  data  evolve  simply  m  time.  Hence,  the  above  RH  and  T  problems 
are  specified  in  terms  of  initial  data;  their  solutions  yield  ?(z;x,y,t) 


and  then  (4.2a)  gives  q(x,y,t). 


We  summarize  the  results  of  (6)  concerning  DSI  (o  ■  1 ,  Proposition 
4. 1.-4. 4)  and  DSII  (c  ■  -i,  Propositions  4. 5-4.0). 

Proposition  4.1  (Bounded  Eigenfunctions)  A  solution  of  (4.2a)  with 
c  *  1  bounded  for  all  complex  values  of  z  *  zR  ♦  iZj  and  tending  to  I 
as  z  •*  •  is  given  by 

Y  +  (  z ;  x .  y ) ,  z:  >  0 

T ( Z ; x , y )  *  /  (4.3) 

f“(z;x,y) ,  Zj  <  0 

where  ti(z;x,y)  satisfy  the  following  integral  equations: 

T*(z:x,y)  »  I  ♦  ^  J*  dC  eiz  ( >' Jj"dnj3n.  <y-n>  IMx-*  >  J  J . 

(»g  +  «t)  (q(£,r|)Y~(z;C,n)  > 

'  r  d{  dnr-  „„  #i-l<y-nii  *  U-OJ1. 

J  X  '  -m  '  -m 

wjfqK.nJY1  <z;C,n>  )  (  4  .  4  ) 1 

(cf.  (3.5)  for  notation). 

♦ 

Assuming  that  the  linear  integral  equations  (4.4)~  have  no  homoge¬ 
neous  solutions,  it  follows  that: 

Proposition  4.2  (Departure  from  Holcrmorphicity )  .  Y  ,  T  are  holomorph 
functions  of  z  for  Zj>0,  Zj <  0,  respectively.  Hence  the  function 

f(z;x,y)  defined  by  (4.3)  is  a  sectionally  holomorphic  function  of  z. 

In  particular,  —  ■  0  for  all  z,  with  zT  ♦  0  and  —  -  Y+(z;x,y)  - 
3  z  1  3  z 

Y  ( z ; x , y )  for  z  *  zR.  The  departure  from  holomorphicity  is  given  by: 

T+(z;x,y)  -  Y  ( z ; x , y )  -  J  dz  1  Y  ( z  * ; x , y  )el Z  +  *  Z  Y 

f (z’ ,z)e'lzJx_lZy,  (4.5) 

for  z  ■  z_,  where  the  scattering  data  f(z',z)  are  given  by: 


f  ,  ,  (2'  ,2) 


■  -  f  dm  f12(2',m)f21<m,z),  f^fz'.z) 

^  -  m 


dn  Qf 22  e 


- ( Z+Z '  ) £+i ( Z“Z '  )  h 


f21(z*,z)  - 


1  f"  ATC  .  .  zm+  -ilzu'lCuU-Dn  - 
17  J  dnX  QTn  e  ,  f. 

—m  *  — m 


•0.  (4.7) 


Proposition  4.3  (Inverse  Problem-Reconstruction  of  the  Potential  q) 
q(x,y)  is  obtained  from: 

q(x,y)  *  -  1  dz'  f  dz¥  (z ' ;x,y)e  f(z  ,z)e  *1 

~~  (4.8) 

where  T~(z;x,y)  solves  the  following  integral  equation: 

,  r-  f  (z-;*,y)eiz'Jxf(z-.*'')«-iz,Jx+i<*"-z,,y 


(  z  ;  x ,  y )  ♦  J  dz"  J  dz  ’ 


z'  -  z  +  iO 


(4.9) 


Finally  from  (4.2b)  we  obtain  the  following: 


Proposition  4.4  (Evolution  of  the  Scattering  Data).  The  evolution  of 
the  scattering  data  from  t  =  0,  f(z',z;0)  is  given  by: 


.  2  2 

f ( z  '  ,  z  ;  t )  -  e~z '  tA30  f(z',z;0)e2  tA30 


(4.10) 


where 


f(z',z;0)  is  given  by  (4.7)  and  A3Q  *  diag(i,-i). 

Proposition  4.5  (Bounded  Eigenfunctions).  A  solution  of  (4.2a)  with 
o  -  - i  bounded  for  all  complex  values  of  z  *  zR  *  iZj  and  tending  to  I 
as  z  -»  •  satisfies  the  following  Fredholm  linear  integral  equation 


¥ ( z  ;  x ,  y )  *  I  +  (G  T ( z ;.,.))( x , y ) 

zr *  * i  *  q 


(4.11) 


where 


i  *  s-Cd'C,I,d’On  '‘"■I?."’ 


(expl  ( m+  i  (  l-J  )z)(*-t)*im(y-n)  )lq(  C.  n)¥(z;  n)  1  >5  y 


(4.12) 


*  MLJ'L  {d-  J_  dn  -  J  'dm/  dn» 

K  1  •  2}  I*-"  “  Km  -m 

(exp(  -  (ro+id+J  .)z)lx-U  +  in<(y-n)  ](qU,n)1'(z;C»n)])^,  <  4  .  I  2  )  2 

ctj  '  wr  c2j  -  :  *  » .2. 

Proposition  4.6  (Departure  from  Holomorph ic i ty )  .  For  every  z  c  C 


0  J  _ 

—  (  z  ;  x ,  y )  *  t(t;x,y)n(z0,zT;x,y) 

3  z  *  1 

where  the  matrix  fi  is  defined  by:  (2^  *  ^21  =  ® 


(4.13) 


ni3  4  j ( z )exp6  —  ( z ; x , y ) ,  i  ♦  j  (4.14) 

Tij(2)  "  17  J  d<J  dnCqtC.nJTIzjC.n))^  expl-ei3(z;c,n)).  i  •  3 

e!2<  z;x,y)  4  2i(xzR  +  yZj),  ^(zjx.y)  *  2  i  ( -xzR+yz  j  )  . 

Propostion  4.7  (Inverse  Problem-Reconstruction  of  q).  q(x,y)  is  ob¬ 
tained  from 


where  Y(z;x,y)  satisfies: 


Finally  from  equation  (4.2b)  we  obtain: 


Proposition  4.8  (Evolution  of  the  Inverse  Data).  The  inverse  data 


at  time  t,  D( zR, Zj ; x,y , t ) ,  is  given  by 


q(x,y)  =  [J,  -jjj-  JJ  r(z;x,y)  n(zR,z^;x,y)dz  *  d  z  ) 

C 


(4.  IS) 


JJf(z;x,y)  f)(zR,zj;x,y) 


(4.16) 


D(  zR ,  Zj  ;  x ,  y ,  t )  ■  exp(z2A3Qt)  D(  zR,  zx  :x,y ,  0  )exp( -z2A3Qt )  (4.17) 

where  Q( zR , Zj ; x , y , 0 )  is  given  by  (4.14)  using  the  initial  condition 
q(x,y,0)  and  Ain  *  diag(i,-i). 


5.  INVERSE  PROBLEMS  IN  0+1 

The  Lax  pair  associated  with  the  PIV  equation 


*  i7(at>2  *  1  y3  +  4ty2  *  2(t2  +  a,y  * 


dt 


is  given  by 


Y  (  2  )  - 

Z 


O  O'! 


0  - 1 


V2 


z  + 


t-(v-e  -e  )  -t  j 

'-y  0  oo  ' 


n 


lo 


( o 


z  + 


u  0  ® 


(5.1) 


'v*  -  ¥  ) 


2  v , 


v37(v“2eo)  -(e(Tv);  J 


Y  (  Z  )  , 


(5.2a) 


Y  (  z  )  . 


(5.2b) 


Indeed  Yzfc  =  Ytj  implies 


dt 

dv 

dt 


=  -  4v  ♦  y2  +  2ty  ♦  4Gq,  ^  =  -u(y+2t), 

2  2  460 

=  '  y  *  9-^  •  y,v  *  (0o  +  e„) y ' 


(5.3) 


where , 

a  -  20^  -  1 ,  8  =  - 80g . 


As  in  the  cases  of  1  +  1  and  2+1,  solving  the  initial  value  problem 
of  PIV  reduces  to  solving  an  inverse  problem  for  Y:  Reconstruct  Y(z;t) 
in  teru  of  appropriate  monodromy  data.  Again  this  inverse  problem  will 
be  solved  in  terms  of  a  RH  problem.  Thus  it  is  essential  to  study  the 
analytic  properties  of  Y  with  respect  to  z.  However,  in  contrast  to 
the  analogous  problem  in  1ST  for  1+1  and  2+1,  this  task  here  is  straight 
forward:  Equation  (5.2a)  is  a  linear  ODE  in  z,  therefore  its  analytic 

structure  is  completely  determined  by  its  singular  points.  In  this 
particular  case  z  •  0  is  a  regular  singular  point  and  z  ■  —  is  an  irregu 
lar  singular  point  of  rank  2.  Complete  information  about  z  ■  •  is 
provided  by  the  monodromy  matrix  and  by  the  Stokes  multipliers 

a,  b,  c,  d.  Solutions  of  (5.2a),  YQ  and  Y)(  normalized  at  zero  and 


S 

L 


O  v  v  ■_  •.  *v 


UK! 


VA".V  Cs'.'.'A'A"  -  .'V  ■.  ,v.'- ."•  A 


infinity  respectively  are  related  via  a  connection  matrix  Eg  with 
entries  ag,  8g,  Yq<  6q.  Taking  into  consideration  the  above  singula¬ 
rities,  there  exists  a  sectionally  holomorphic  function  Y,  with  jumps 
across  the  four  rays,  arg  z  ■  -  j,  j,  and  with  singularities 

at  z  »  0,  z  ■  - .  The  jumps  are  specified  by  a,  b,  c,  d  and  the  nature 
cf  singularities  by  Mg,  M^.  This  leads  to  a  matrix,  singular,  dis¬ 
continuous  RH  problem,  defined  on  the  above  rays  and  specified  in  terms 
of  the  monodromy  data 


Monodromy  Data  (MD) 


{ a , b, c, d ,  aQ  « 8g  <  Y0  *  ^ • 


A  consistency  condition  of  the  above  RH  problem  yields 

<  i  VM.  *  EoVEo' 

where  are  the  Stokes  matrices  uniquely  defined  in  terms  of  the  Stokes 
multipliers.  Using  (5.5)  and  certain  similarity  arguments  it  can  be 
shown  that  all  MD  can  be  expressed  in  terms  of  two  of  them.  Furthermore, 
equation  (5.2b)  implies  that  the  MD  are  time  invariant.  Hence  the 
above  basic  RH  is  specified  in  terms  of  two  initial  parameters  (these 
two  initial  parameters  are  obtained  from  the  two  initial  data  of  PIV). 
The  solution  of  this  RH  problem  yields  Y(z;t)  and  hence  (5.2a)  yields 

y  ( t )  . 

The  above  basic  RH  problem  can  be  simplified  considerably:  (l) 

Assume  0  0Q  <  1,  0  <_  8^  <  1  •  8g  *  2'  tben  the  above  RH  problem  is 
regular.  It  is  interesting  that  the  basic  RH  problem  can  be  used  to 
obtain  Schlessinger  transformations  which  shift  eQ  and  ©^  by  a  half¬ 
integer.  By  using  these  transformations  the  general  case  is  reduced 
to  the  regular  case.  (ii)  The  basic  RH  problem  can  be  mapped  to  a 
sequence  af  two  RH  problems,  one  on  the  line  arg  z  **  ^  and  the  other 
on  the  line  arg  z  *  -  i.  The  first  one  is  continuous  (both  at  x  =  0 
and  x  ■>  «•)  ;  furthermore,  it  can  be  solved  in  closed  form.  The  second 
one  is  discontinuous  both  at  x  *  0  and  x  =  By  using  standard  auxi¬ 

liary  functions  one  maps  the  discontinuous  problem  to  a  continuous 
one.  Then  the  theory  of  continuous  RH  problems  on  simple  contours 
can  be  used  to  establish  uniqueness  and  existence  of  solutions.  Elemen¬ 
tary  solutions  of  PIV,  expressible  in  terms  of  Weber-Hermite  functions 
are  natually  obtained  within  the  above  formalism.  w@  summarize  the 
results  of  (7)  concerning  PIV. 


Q { z )  *  Diag(q.-q),  q(z,t)  4  ♦  zt,  Y„(z)  is  holomorphic  at 


(iii)  Yq ( ze  " )  =  Yq(z)M0,  Mq  * 


xp ( 2 1 )  2  in Jexp( 2 ih0q ) 


exp( -2in6  ) 


,(5.4) 


7m  0  if  V  I<  J"‘f  ®Q  *  i- 


(iv)  Y  2 ( 2 )  «  Y  ^ ( z  )G  j ,  Y  3 ( z  )  -  Y2(z)G2,  Y4<z)  =  YjtzIG^, 

Y  ^  (  z  )  «=  Y4  (  ze21"  )G4M-( 


(5.5) 


where 


.  g2  * 


a  1 


1  d 


0  1 


0  1 


.  G3  4 


1  0 


c  1 


,  M  4  exp( 2 l nD  ) . 


(5.6) 


(Iv)  Y  j  (  z  )  -  Yq  (  z  )  Eq  ,  Eq  * 


a0  B0 


Y0  40 


,  det  Eq  =  1 


(5.7) 


Furthermore,  the  parameters 


MD  4  (a,b,c,d,  aQ ,  Bq,  Yq.  6 Q > 


satisfy  the  following  consistency  condition. 


(vi)  (  n  G  . )  M 
j-)  3  * 


Eo'Mo'Eo- 


(5.9) 


Proposition  5.2  (Properties  of  Monodromy  Data ) 

(i)  The  monodromy  data,  MD,  given  by  (5.8)  and  defined  in  Proposition 
5.1,  are  time- invr lant . 


,TiPi2in 


11 )  All  of  Che  MD  can  be  expressed  in  terms  of  two  of  them. 


(iii)  ( 1 *bc)exp( 2ix0 )  +  (ad  ♦  ( 1 +cd ) ( 1 ♦ab ) ]exp( -2 i w0  )  ■  2  cos2n0. 

•  •  *D 

(5.10) 

In  what  follows  we  formulate  a  RH  problem  for  the  case  that 
0  <  0q  <  1/  0  <  0^  <  1 .  This  assumption  leads  to  a  regular  RH  problem. 
The  general  case  follows  by  considering  this  result  and  Schlessinger 
transformations. 


Theorem  5.1  (Inverse  Problem).  Consider  the  following  matrix,  regular 
homogeneous  RH  problem  along  the  four  rays  Cj,...,C4  (Figure  5.1): 


Determine  the  sectionally  holomorphic  function  T(z),  T(z)  =  TMz)  if 
z  is  in  Sj,  j  =  1 , .  .  . ,  4 ,  from  the  following  conditions: 


1.  f j  satisfy  the  jump  conditions 


'*2{c-)  m  ^  t  ( C  <  C  >  -  *  3  ( C  )  *  *2(c,g2(c)'  iVc)  *  *3(C>93<C> 

2in 


*,<C>  =  Vc  •  >94<C> 

along  the  rays  ,  Cj,  C4 ,  Cj  respectively,  where 


(5.11 


91  *  e°G3e"Q- 


1,2,3,  g.  4  eQG.e"QM  . 

4  4  ® 


2.  f  <z>  -  (-)  “(I  +  0 ( — ) )  as  I z I  -  -. 

z  z  1  1 


(5.12) 

(5.13) 


3.  f(z)  has  at  most  an  integrable  singularity  at  the  origin  with  a 
munodromy  matrix  given  by 


T,<ze21")  -  »,I*Ie;V0-  2 


0  . 


(5.14) 


In  the  above,  G^,  Q,  M^,  D^,  Mg  are  defined  in  Proposition  5. 


4.  The  monodromy  data  MD,  given  by  (5.8),  satisfy  the  properties  given 
in  Proposition  5.2(ii).  Then: 

(i)  The  above  RH  problem  is  discontinuous  both  at  the  origin  and  at 
infinity.  Actually 


- 1  - 1 

-  Eq  Mq  E0.  Z  -  0; 


(1  g  -  M  ,  z  -►  ». 

’  h  am 


(5.15) 


r 


3.  Klz)  has  at  most  an  integrable  singularity  at  the  origin  with  a 
monodromy  matrix  given  by 


K(ze2l1,>  -  K(z)h]  (OJEQ’MQEQh"1  (0)  ,  z  -0. 


(5.18) 


The  above  RH  is  discontinuous  both  at  the  origin  and  at  ir.cinity. 
Actually  if  gR  ,gR  denote  the  jump  matrices  along  Cj,  respectively 
then  1  3 


g^gK  -  hj  ( 0 ) Eq  'hq  1  EQh j 1  ( o ) ,  z  0;  g"^  ~  *m.  z 


(5.19) 


However,  the  above  RH  problem  can  be  mapped  to  a  continuous  one  using 
the  auxiliary  functions 

±e_ 


'ztl  '  '  'zil 


to  remove  the  above  singularities. 


¥  is  related  to  K  via: 


(5.20) 


¥  =»  Kh  if  z  in  Sj+Sj;  ¥  =  KhM,  M  v  Diag(l,-a/c) 


(5.21  ) 


if  z  in  S^+S^ 

(i.e.,  ¥ |  ■  Kth1'  T2  *  K1h2r  *3  =  K2hiM'  t4  =  K2h1M)‘ 

Proposition  5.3  (The  Solution  of  PIV).  Let  ¥(z)  be  the  solution 
matrix  of  the  inverse  problem  formulated  in  Theorem  5.1.  Then  y(t). 


y ( t )  -  -(^  gr  ♦  2 1 ) ,  u  *  -21im  ¥, , ( z  )e‘2q< z  ’  , 

1*1  — 

solves  PIV. 


(5.22) 


6.  INVERSE  PROBLEMS  IN  n  SPATIAL  DIMENSIONS,  n>  2 

Consider  the  inverse  problem  associated  with  the  following  system  of 
N  first-order  equations  in  n+1  dimensions: 


¥  +o  t  J.¥  -  q¥,  o  *  oB  +  ioT,  o.  ♦  0,  n  >  1, 


X0  t-1  1  Xl 


(6.1) 


9 

i 

& 


where  q(Xg,x)  is  an  N  *  N  matrix-valued  off-diagonal  function  in  Fn* 
decaying  suitably  fast  for  large  xQ,  x,  and  the  are  constant  real 

diagonal  N  *  N  matrices  (we  denote  the  diagonal  entries  of  J,  by  jl, 

N  ^  ^ 

.  ..,J^).  Alternatively,  using  the  transformation 


f(z,xQ,x)  -  u(z,xQ,x)exp( i  I  z t -  oxqJj)],  z  cCn, 


(6.2) 


equation  (2.13)  becomes 


a  l  (J.WV  ♦  is.  Ut,wJ ) 

t«i  1  *  1  * 


(6.3) 


We  assume  that  n  ^  N,  otherwise  the  entries  of  the  matrices  will  be 
linearly  related  and  one  can  always  reduce  n  by  a  change  of  coordinates. 
An  inverse  problem  in  this  case  is  defined  as  follows:  Given  appropriate 
inverse  data  T,  where  T  is  an  N  x  N  matrix-valued  off-diagonal  function 
of  suitable  inverse  parameters-  reconstruct  the  potential  q.  Again 
there  exists  a  u  which  is  bounded  for  all  complex  z,  z  cCn.  3u/3z 
depends  on  appropriate  inverse  data  T( z„ , z . ,m- , . . . ,m  ) ,  z„eFn,  zTeF? 

im  J*,  K  I  &  n  K  1 

m  c  R.  T  satisfies  -il — =  , *Lg.  .Using  this  equation  and  introducing 


Born  variables,  *Zj 


Z  ,  W  CO 


—  — ,  0  efl  ”  1 

,,w,x;  Wq  £  F  ,  w  c  F  ,  \  e  C 


(6.4  ) 


one  obtains  a  characterization  equation  for  the  inverse  data: 

dyT>  dXP  NjjlTl(w0.w,XP'  ) 

Tlj(w0,w)  *  Tlj(wQ,w,  X)  -  \  \  2  - - - 1 - 

*  Xp  - 


,  (6.5) 


where  N  is  a  quadratic  function  of  T.  That  is,  Ti;3(z,m)  is  appropriate 
inverse  data  iff  the  right-hand  side  of  (6.5)  is  independent  of  X. 

Hence,  equation  (6.5)  serves  as  both  characterizing  T1*  and  defining 
T*^.  This  equation  was  first  introduced  by  Nachman  and  Ablowitz  ( 8 ] . 
Using  equation  (6.5)  and  taking  the  limit  of  u  as  Ixl  ■*  •  we  show  that 
the  general  problem  of  reconstructing  an  N  x  N  potential  q  in  n+1  spatial 
dimensions,  is  reduced  to  one  of  reconstructing  a  2  «  2  potential  with 
entries  q*^,  qJ*  in  two  dimensions.  The  inverse  data  needed  for  this 

»  >  i  •  •  , 

reconstruction  is  precisely  T1  ,T^ 1 .  This  reduction  makes  crucial  use 
of  the  existence  of  redundant  scattering  parameters.  In  this  sense  it 
is  the  analog  of  the  Born  approximation.  However,  the  crucial  difference 


is  that  while  in  the  inverse  scattering  of  the  multidimensional  Schro- 
dinger  equation  one  can  reconstruct  the  potential  in  closed  form,  here 
one  can  only  reduce  the  general  problem  to  one  for  2«  2  matrices  in 
two  dimensions.  This  reduced  problem  was  solved  In  (6J.  In  the  follow¬ 
ing,  we  summarize  the  results  of  (9). 

Proposition  6.1  (Bounded  Eigenfunctions).  The  f unction  yi  ( xQ , x , z ) 
defined  below,  solves  equation  (6.3),  is  bounded  for  all  complex  values 
of  z  and  tends  to  I  for  large  k: 


u 1 ^ ( Xq , x , z  )  «  4  *  ^  ♦ 


sgn  ( c  j  J , ) 


1  r  exp( ifi 


exp( i0 1 ] (Xq-Cq/X, ,  ,z) ) 


)  -  oj;(x0-c0) 


lqu)i ,x2-(x1-51)jJ/j| . xn*,x,-51 cC"' ,6-5) 

where  6 1 ^  is  defined  by 


0  ;i(x0.x1,z)  * 


xo  !°  I  ziT 


xi(ozt,i 


,z»=  7 i  *lzi 
1  K  lI 


(6.6) 


Equivalently  satisfies 


v  J(xn,x, z)  *  6 


‘j  ,  sqnlqIJ  1 1  _  f  d(  dt|c  (■ 

J„n+1  -0  s  n-U_n-i 


exp(  ifli^(x0-^0,ic1-t1,z)  KquJ^Kg.C.z)  , 
X1 1 -OJ1 *  x0~^0  * 


(6.7) 


where 


i  n 

♦  dm-...dm  .  a  (x,m)  ♦  E  m.fx.-x,  — <■) ,  c  ♦ 
2  n  t.2  1  1  1  j‘  n 


(6.8) 


Proposition  6.2  (Departure  from  Holomorphicity ) .  Let  yt  J  he  defined 
by  eq .  (6.5).  Then 


(6.9) 


^  (x.,x,z)  ■  I  Y1(J^  *  J^)exp(  16  1 3  (  xQ  ,  X  1  ,  z  )  ] 

a*p  i'3 

«  c.  ,  f  „  .  dm2exp{ ia * ( x ,m) IT1 ^ ( z ,m) o( x. , x , A* ^ ( z ,m) )E. , , 
n-i  JRn-i  oi] 

where  S^(Xq,Xj,z),  a*(x,m)  are  defined  by  (6.6),  6-8)  respectively;  Et 


is  an  N  «  N  matrix  with  zeros  in  all  its  entries  except  the  i]th,  which 
equals  1;  and  X  ^  and  T*^  are  given  by 

,i 


J 


n  J 


x|^(z,m)  $  ,  z .  ),  X^tz.m)  =  (z  +m  ,  z  );  r*2,..,n 

R  ‘“2  j|  ’i  r  rR  r) 


Y  1  *  o/4*  i|  jJox|  , 

Tlj(z,m)  *  f  n+1  dCQd^  expJ-iB1  J(  CQ,  C,  ,z)-ial  (  C,m)  ]  (qu)13(  CQ,  C.z)  • 

J  K 


(6.10) 


Proposition  6.3  (Characterization  of  T) 


(a)  Assume  that  3u/3zp  is  given  by  Eq .  (6.9)  and  the  T1-,(z,m)  is  given 
by  (6.10).  Then 


L*pT* ^  (  z  ,  m  )  = 


n  .  2 

l  c  .  .  dM  T 

t-i  J*"-' 


( j1-j j) ( j1-jl  ) 
r  r  p  p 


ll(Xl  ^(  z,M)  ,m-M)T  t;,(z,m) 


*  pi T 1 ( * . m ) , 


(6.11) 


where 


l*2  *  <Ji  -  j2)  - - (j2  -  J?)  - 


rp 


(6.12) 


32. 


3  z 


(b)  Assume  that  3u/32  is  given  by  Eq .  (6.9)  and  that  3  u/3z  3z 

p  r  p 

is  symmetric  with  respect  to  r,  p.  Then  Tl^(z,m)  solves  (6.11). 


Following  A.  Nachman  and  M.  J.  Ablowitz  we  introduce  appropriate  j 

Born  variables.  Then  equation  (6.11)  can  be  integrated.  Furthermore,  s 

we  can  compute  the  limit  of  T1^  in  the  new  coordinates  as  |xp|  •*  •  j 

( see  below) :  } 

? 

i 

i 

\ 

> 

i 

i 


Let  Wq  j ,  wJJ,  wfc/  t  =  2 ,  .  .  .  » n  e  *  and  X|lcC  ,  l  *  2 . n,  be 

defined  by 

n  J1-^  ^  4  n  n  J1 

w*3  4  £  r-  |  o  |  2  z  ,  w}3  +  -  £  —  ,j£-(oz  )  -  £  m  -4, 

0  r.l  «I  rI  '  r-1  CjjJ  r  1  r-2  c  j| 


♦  mt# 


-,  t  »  2 ,  .  .  .  ,  n . 


(6.13) 


Assume  that 


(J^-J3) ( J^-J3 )  ♦  (  J*-J?  )  ( Jr-J3) ,  for  all  distinct  i,j,r  and  p*  1 . 

1  1  p  p  1  1  p  p 


(6.14) 


For  convenience  of  writing  we  usually  suppress  the  superscripts,  i,j 

in  Wq,w1#  x-  Let  *  denote  m  denote  m2,...,mn,  x  denote 

y-,,...,  y  ,  w  denote  w,,...,w  .  Then  we  have  the  following. 

A  i  A  n  i  n 

(a)  The  inverse  of  the  transformation  Z,m  -  «Q .  w»  x  is  given  by 


! t  ”xt(Ji-Jl)'  mt  “  V  1  *  2 . .  2i  m  + 


(o/|o|2)wq  wgJ* 

j  3  _  j  i 

J,  J, 


(6.15) 


(b)  In  the  new  coordinates,  Eq.  (6.11)  with  r  =  1  becomes 

~ — ( w.  ,  w ,  x  )  -  N*3(T)  (w.  ,w,  x  )  ,  P  *  2  ,  .  .  .  ,  n  .  (6.16) 


(c)  In  the  new  coordinates, 

T 1 ^ ( wq  #  w , x )  *  |  n+1  exp( -i (w0C0+wC ) I (qu ) 1 ] (t0,t,w0,w,x) , 

where 


£  w  £ 

r-1  r  r 


(6.17) 


V-V  '  v  ^ ^ 


yri.Pt.n^JHLWJ "i J 


(d)  Let 


gXj  ♦  Ul:j(x0/x.w‘J,wiJ,  y/J),  fir 


O  ij  vij; 


Then 


the  gX^  satisfy 


lim 

IXph- 


fi^(x0,x,wQ,w) 


agnfOjjj ) 


2s  i 


'n-1  Jp2n 


dx^dx'dw  exp(  i  { (xQ-x^)w0  -Mx-x’lw}) 


(  Xq  ,  x  *  >  0 1  ■*  (  *q  > x '  ,  wQ ,  w )  * 


xrx;-  aj\[xQ-  X*) 


^ i ^ ( x0 ' X ' Wg  *  w  ^ 


*  1  ♦ 


sqnCOjJ*  ) 

2s  i 


n-1  j_ 


2n 


dx^dx '  dw^q  1  (  x^ ,  x '  )  g  { x£ ,  x '  ,  wft ,  w) 


x,-x;  -  oj3(x0-x') 


— —  ,  yX^  *  0,  for  all  l, 

(6.19) 


(e)  lim  T1  ^  ( wQ , w , x )  *  J  n+,  d^d^  exp[ - i ( wQC0+wt ) ]  * 

lxph-  * 

q1  j(C0.OwJj(C0^.w0,w)  ♦  Tlj(wQ,w).  (6.20) 

(f)  The  basic  characterization  equation  is  given  by 


T1 ^ ( wQ , w ) 

m  Tlj(w  ,w,X)  -  ~  f  2 

dXp  dXp  n)^(T)(w0,w,xP  ) 

®  I  ■  r  6 

x  -  x 1 
*p  *p 

~  Jv. 

where  x*3 

denotes  X2**’-»  Xp_,» 

XP'  Xp+1 .  V 

It  follows 

from  the  above  that  as 

|x_|-»«>  the  yl^’s  decouple. 

Furthermore,  the  y 1  ^ ,  y?^  satisfy  a  system  of  two  equations  depending 

i  1  i  i  1  1 

on  q  J,  qJ  .  It  turns  out  that:  (a)  By  introducing  appropriate  spatial 
variables  the  y|^,  yj^satisfy  equations  in  two  spatial  £invttiiions . 


.  s  .if.  <.  *  - 

'a  /  V'x*.  /.V. /.V  '/ 


(b)  The  inverse  data  needed  to  reconstruct  u  x  >  u L  (and  hence  q 
q^)  can  be  obtained  from  T13: 


Proposition  6.4  (Reconstruction  of  q).  Let 


J2JJ  ‘  JrJ2 

a,  *  -M - S, 


JtJr  - 


z  |if3  _  7D7*  z  —  t  D  j  i 

J1J2  J1J2  J1J2  J1J2 


,  r  *  1  ,  .  .  .  ,  n , 


(6.22) 


where  for  convenience  of  writing  we  have  suppressed  the  dependence  of 
a r  /B  r  on  i ,  j .  Let  (q  t  F,  5  e  *n. 


x0  *  C0>  xi  =  ^ i '  x2  =  ^ 2 ' 


6.23) 


Xt  =  Cl  *  atC1  +  8t  *2'  1  *  3 . n- 


Then  we  have  the  following: 


(a)  The  system  (6.19)  becomes 


■  =9"  ^  /  J  asjdtju.-tl-  oJt(5„  -  CJII 


J  ^ 

exP(  idlju0-5o*5  r^i  >i]  iqA  i  ( ^6  ^  i  '^2  '  i  i  ’“T'^3 . vz 


03j(Co.C.i)  «  ,  >  sgn  -4^  J  d^dc;(crt;  -  ojJUq-^)]-1 


-  ^  3 . V£)' 

where  ' 

n  J  j  “  J  i  .  ,  , 

z  4  I  (z  a  ♦  — - - j-  z  6  ),  8  J(x.,x.  ,zl 

r-1  J?  -  J, 

.  .  J  ?  ,  (a z) 

813<x0#x1,z)  *  — - - 1  xo  I  a  I  zx  “  *,  - —  I- 

1  3  . 


(6.24) 


(6.25) 


(b)  T  J  in  the  new  coordinates  becomes 

Tlj(z,m)  »  j  n+)  dZ'dZ1  exp( - 1 8 1 3 ( Cq , Z ;  .  z  ) 


(6.26) 


*  'V')  '  M  “1 1  *  £, 


) , 


r-3 


where 


m,  ♦  m,  ♦  £  mg,  m  -  m. 

i  i  r_3  r  r  l 


l  -  3  ..... n . 


(6.27) 


< c )  The  inverse  data  associated  with  (6.24)  and  the  analogous  problem 
for  Uj*.  Uj*  »re  given  by  T*  ^ ,  f ^ .  Let 


J"  j  * 

'r‘3(2-£2-'t  TT~*  «3 . <„>  *  c„-,  L-,  ..pliSjKj-t,  -f> 

J1  F 


♦  i  I  m  c  |T13(£,m) 
r-3  r  r 


(6.28) 


Then 


-|,<j . t„i  ♦  J  2  expi-iSi:luj.t;.jii 


iq13C3]>(c;.ij.t2  -  <«,-«;>  -f.  h . 

J1 


(6.29  ) 


Equations  (6.1)-(6.3)  with  a  -  -1  lead  to  a  system  which  appears 
to  be  physically  more  interesting:  (a)  Since  the  system  is  hyperbolic  one 
may  consider  the  physically  important  question  of  inverse  scattering  (IS); 
i.e.,  given  a  scattering  amplitude  function  S(X,k)  find  the  potential 
q(Xg,x).  (b)  A  special  case  of  the  above  system,  namely  if  the  J^s  are 
constrained  by 


Jn  - 

-? — e_ 


_e - E 


r  *  Jr 


Jr  -  Jr 


p.  r  -  1 .....  n ,  1 , j. t  =  1 . N,  (6.6) 


is  associated  with  the  N-wave  interaction  in  n+1  spatial  and  one  temporal 
dimensions  [10].  The  above  system  can  be  considered  as  a  limiting  case 
of  16. 1 )  —  (6.3)  [8].  Alternatively,  it  can  be  considered  on  its  own  right 
(11);  the  problem  of  reconstruction  can  be  reduced  to  one  for  a  2  *  2 
matrix  problem  in  two  spatial  dimensions. 


ACKNOWLEDGEMENT a 

This  work  was  partially  supported  by  the  Office  of  Naval  Research  under 


Grant  Number  NOOO 14-76-C-0867  and  the  National  Science  Foundation  under 
Grant  Number  DMS-8501325.  A.  S.  Fokas  thanks  Professor  Lakshmanan  and 
his  school  for  their  warm  hospitality. 

REFERENCES 

[1)  R.  G.  Newton,  Scattering  Theory  of  Waves  and  Particles  (Springer, 
Berlin,  1 982 ) . 

( 2 1  F.  0.  Gakhov,  Boundary  Value  Problems  (Pergamon,  (966);  N.  P.  Vekua, 
Systems  of  Singular  Integral  Equations  (Gordon  and  Breach,  1967). 

(3)  S.  Prossdorf,  Some  Classes  of  Singular  Equations  ( North-Holland , 
1978) . 

[4]  L.  Hormander,  Coaiplex  Analysis  in  Several  Variables  (Van  Nostrand, 
1966). 

(51  R.  Beals  and  R.  R.  Coifman,  Comm.  Pure  Appl.  Math.  37  (1984)  39. 

(61  A.  S.  Fokas,  Phys.  Rev.  Lett.  51  (1983);  A.  S.  Fokas  and  M.  J. 

Ablowitz,  J.  Math.  Phys.  25  (1984)  2505. 

(7)  A.  S.  Fokas,  U.  Mugan,  M.  J.  Ablowitz,  A  Method  of  Solution  for 
Painleve  Equations:  PIV,  PV,  preprint  INSI73  (1987). 

(81  A.  I.  Nachman,  M.  J.  Ablowitz,  Stud.  Appl.  Math.  71  (1984)  251. 

(91  A.  S.  Fokas,  J.  Math.  Phys.  27  (1986)  1737. 

(101  M .  J.  Ablowitz  and  R.  Haberman,  Phys.  Rev.  Lett.  35  (1975)  1125. 

(11 j  A.  S.  Fokas,  Phys.  Rev.  Lett.  57  (1986)  159. 

(12)  A.  S.  Fokas  and  M.  J.  Ablowitz,  Comm.  Math.  Phys.  91  (1983)  381. 


a.v.  M;khaiUv  3D(»<?8i)  ia,i(i3-in 


WWW 


■\r. 


•i  . 


*  -  «.*.  »*. . 


lA WjJ 


